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Simulations with Complex Measures 
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Towards a solution to the sign problem in the simulations of systems having indefinite or complex-valued 
measures, we propose a new approach which yields statistical errors smaller than the crude Monte Carlo using 
absolute values of the original measures. The ID complex-coupling Ising model is employed as an illustration. 



1. The sign problem 

The simulation of many interesting statistical 
and quantum systems suffers from the sign prob- 
lem when the generating function measures are 
indefinite or are complex [1]. Some important ex- 
amples are Lattice QCD with finite temperature 
and density and chiral gauge theory. 

If the measure p(x) of a generating function suf- 
fers from the sign fiuctuation then another posi- 
tive definite function p(x) must be chosen for the 
Monte Carlo (MC) evaluation of the expectation 
value of an obsevable 0, 

(6) = j{Q{x)p{x)/p{x))p{x) I {{!)), (1) 
((1)) ^ [(p(x)/p(x))p(x), (2) 

In general it is desirable to choose p indepen- 
dently of 0, and we will concentrate on the 
estimate of the denominator ((1)) because of 
its appearance in all the measurements. It 
is a simple variational problem to show that 
the MC probability density, properly normalised, 
which minimises the variance of ((1)) which is 
(T/-y/:#independent configurations and where 

^2 = / \p{x)lp{x)-{{l))tp{x), (3) 

J X 

must be 

p{x) = \p{x)\ I jjp{x)\. (4) 

Thus the sign of p{x) is now treated as part of 
the quantity whose expectation is to be measured; 
hence the name average sign approach [1]. 



However, when the denominator ((1)) is vanish- 
ingly small, cr^, though minimised, is 1 ^ ((1)) 
since pl\p\ = ±1. Then the evaluation of (0) be- 
comes unreliable unless the number of indepen- 
dent configurations is many orders of magnitude 
greater than the large number The fiuc- 

tuation of sign of the measure over configuration 
space thus renders ineffective the sampling guided 
by this crude MC method (4). This is the content 
of the sign problem. 

Many approaches have been proposed to tackle 
the sign problem but none is satisfactory [2]. In 
the next section we will improve on the one pre- 
sented above. 

2. A new approach 

We can write the generating function as 
integrals over two configurational subspaces, 

p{x) = Iy p{^j ^)) ill such a way that the 
multi-dimensional integral over Y can be evalu- 
ated analytically or well approximated: 

eix) = J^p{X,Y)- (5) 

As in the last section, one can easily prove that 
the MC weight g that minimises 

a'-' = [ \e{x)re{x)-{{i))t~e{x), (6) 

Jx 

is 

~e{x) = \eix)\ / j^\e{x)\ . (7) 

It then follows that the variance for this new 
weight is not bigger than that for p{X,Y). 
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the second line follows from (4) and (7); the last 
line is from (5) and always less than or equal 
to zero because of the triangle inequality. The 
equality occurs if and only if p{X, Y) is semi- 
definite (either positive or negative) for each X- 
configuration. In particular, when there is no sign 
problem in the first place, expression (6) yields 
the same statistical deviation as the crude one. 

To deal with complex integrands, of which in- 
definite measures are special cases, we adopt the 
definition (3) of the variance extended to the ab- 
solute values of complex numbers. Statistical 
analysis from this definition is the same as the 
standard analysis; except that the range of un- 
certainty should now be depicted as the radius 
of an 'uncertainty circle' centred on some central 
value in the complex plane. All the manipulations 
above remain valid. 

Our approach is now clear. The measures are 
first summed over a certain subspace, the inte- 
gration (5) above, to facilitate some partial phase 
cancellation. Absolute values of these sums are 
then employed as the MC sampling weights (7). 
The numerator in (1) can be obtained from an 
appropriate derivative of the generating function 
after the partial summation. 

3. Complex coupling Ising spins 

Owing to the nearest-neighbour interactions, 
the lattice can be partitioned into odd and even 
sublattices, of which the Ising spins Si (= ±1) on 
site i(G the sublattice) do not interact with each 
other. Absolute values of sums of the complex- 
valued weights over the even sublattice, say, are 
the new Monte Carlo weights. 

In our simulations, periodic boundary condi- 
tion is imposed on the ID chains of 128 Ising spins 
which become 64 spins after the partial summa- 
tion. Ensemble averages are taken over 1000 con- 
figurations, out of 2^^® possible configurations. 



They are separated by 30 heat bath sweeps which 
is sufficient for thermalisation and decorrelation 
in all cases except perhaps one, as will be demon- 
strated shortly. A heat bath sweep is defined to 
be one run over the chain, covering each spin in 
turn. The numbers of trials per sweep are dif- 
ferent before and after the partial summation be- 
cause the numbers of spins to be updated are not 
the same. Both hot and cold initial configurations 
are used and will be indicated when needed. The 
crude MC weights of the average sign method are 

Pi{s}) ~ n exp(JreSiSi + l) (9) 

all sites 

and the improved MC weights are 

Qi{s}) ~ n |cosh(J(s,+s,+i))|, (10) 
odd sites 

with coupling J = { J re, Jiva) and no external field. 

We have presented some measurements in the 
table by exact, improved MC and crude MC 
methods respectively. Expressions for the mag- 
netisation and susceptibility are obtained from 
appropriate derivatives of corresponding partition 
functions. The complex correlation length ^ is 
calculated from the transfer matrix eigenvalues. 

The autocorrelation of the unit operator 
in Fig. 1 show that the noise is too overwhelm- 
ing in the crude simulation to tell whether 30 
sweeps are sufficient for thermalisation or not. In 
contrast, the improved simulation is very well- 
behaved. 

4. Conclusions 

Owing to the sign cancellation in the partial 
sums, our approach can offer substantial improve- 
ments over the crude average sign method and 
may work even when the later fails, in the region 
of long correlation length and vanishing partition 
function. 

The choice for splitting of the integration do- 
main is arbitrary and its effectiveness depends 
on the physics of the problem. If the inter- 
actions are short-range (not necessary nearest- 
neighbour), maximal, non-interacting sublattices 
can always be chosen to provide a natural split- 
ting, which is the particular splitting for our il- 
lustrative example. The approach of Mak in [2] 
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Coupling 


Magnetisation per spin 


Susceptibility per spin 




(0.1,0.1) exact 
improved 
crude 


(0,0) 

(0.0007,0.0000) [0.0024] 
(-0.0032,0.0016) [0.0062] 


(1.1971,0.2427) 
(1.1921,0.2459) [0.0356] 
(1.1651,0.2883) [0.1460] 


(1.96,0.78) 


(0.01,0.1) exact 
improved 
crude 


(0,0) 

(-0.0006,-0.0001) [0.0021] 
(0.0023,-0.0080) [0.0055] 


(0.9999,0.2027) 
(1.0481,0.2231) [0.0268] 
(1.0868,0.1374) [0.1198] 


(2.29,1.47) 


(0.01,0.5) exact 
improved 
crude 


(0,0) 

(-0.0031,-0.0047) [0.0027] 
(-0.0522,0.0820) [0.2074] 
(0.0123,-0.0723) [0.0830] 


(0.5512,0.8585) 
(0.5604,0.8163) [0.0428] 
(0.6808,5.4507) [9.2788] 
(1.5323,1.3885) [2.3122] 


(0.60,1.55) 




Figure 1. The autocorrelation versus the number of sweeps for J = (0.01,0.5) from a total of 2 x 10® 
configurations. The boxes are from crude weights; triangles, improved weights. 



could be considered as yet another particular and 
non-trivial splitting where the inner integration 
was approximated in a certain manner. 
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